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ABSTRACT 

We calculate the volumetric rate of superluminous supernovae (SLSNe) based on 5 
events discovered with the ROTSE-IIIb telescope. We gather light curves of 19 events 
from the literature and our own unpublished data and employ crude k-corrections 
to constrain the pseudo-absolute magnitude distributions in the rest frame ROTSE- 
IIIb (unfiltered) band pass for both the hydrogen poor (SLSN-I) and hydrogen rich 
(SLSN-II) populations. We find that the peak magnitudes of the available SLSN-I 
are narrowly distributed (M = —21.7 ± 0.4) in our unfiltered band pass and may 
suggest an even tighter intrinsic distribution when the effects of dust arc considered, 
although the sample may be skewed by selection and publication biases. The presence 
of Oil features near maximum light may uniquely signal a high luminosity event, 
and we suggest further observational and theoretical work is warranted to assess the 
possible utility of such SN 2005ap-like SLSN-I as distance indicators. Using the pseudo- 
absolute magnitude distributions derived from the light curve sample, we measure the 
SLSN-I rate to be about (32+ 26 ) events Gpc -3 yr _1 at a weighted redshift of ~Z = 
0.17, and the SLSN-II rate to be about (lBltH 1 ) events Gpc -3 yr _1 /if 1 at z = 0.15. 
Given that the exact nature and limits of these populations are still unknown, we 
discuss how it may be difficult to distinguish these rare SLSNe from other transient 
phenomena such as AGN activity and tidal disruption events even when multi-band 
photometry, spectroscopy, or even high resolution imaging are available. Including 
one spectroscopically peculiar event, we determine a total rate for SLSN-like events of 
(199±86 7 ) events Gpc" 3 yr" 1 at z = 0.16. 

Key words: supernovae: general. 



1 INTRODUCTION 

There is a growing sample of supernovae with peak lu- 
minosities over 30 times brighter than average (based on 
M = —16.8 fro m the volume limited Lick Observatory Su- 
pernova Search, iLi et alJlioTlTl . These superluminous super- 
novae (SLSNe), were not identified in the first 60 years of 
studies following the pioneering work of Zwicky and Baade, 
which marked the be ginning of sys tematic searches for ex- 
tragalactic transients (|ZwickvHl938T ). 

This lack in what would naively seem to be the eas- 
iest kind of supernova to discover could in part be ex- 
plained if the historical searches were simply looking in the 
wrong place; indeed, the sample of SLSNe publis hed thus far 
show s a preference for low luminosity galaxies (|Neill et al.l 
120111 ') and possibly galaxy cores-two environments neglected 
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by early surveys. But the greater distance over which the 
SLSNe are visible further suggests that the intrinsic rates 
must be low lest they be discovered in the background of a 
targeted galaxy. 

The origin of these events is still a matter of debate 
(for a recent review, see lGal-YarrJ|2012h . While the host 
environments and energetics may best be understood in 
the context of massive stellar explosions from young, star- 
forming environments, the ultimate power source is less 
clear. The decay of radioactive 56 Ni, the source of a Type 
la explosion's brillian ce, may be i nvolved in certain e vents 
jGal-Yam et al.ll2009l . but see also lDessart et al.ll2012T ), but 
for others, it cannot b e more than a minor c ontributor (e.g. 
Pastorello et all |2010| ; lOuimbv et~ai1 l201ll ; IChomiuk et al.1 
20111 ). Some events show narrow emission lines that indi- 
cate that an interaction with a slow-moving CSM may drain 
the large store of kinetic energy from th e supernova ejecta 
and transfer this into radiated light (e.g. ISmith et al.ll2007l ; 
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lOfek et all 120071 ; ISmith et all l2008al ; brake et al.l l2010h . 

Othe rs show no such e v idence for a long-l ived interaction 
(e.g. IMiller et all 120091 ; iGezari et all 120091 ). Some SLSNe 
may draw their power from a compact object that forms 
in the collapse of the progenitor star (e .g. iKasen fc Bildsten] 
l20ld : IWooslevll2010l : IOuved et alJl2012h . 

Whatever their origin, the high luminosities of SLSNe 
make them detectable from large distances, and they are 
thus capable of shepherding information from the early uni- 
verse to the present. Intervening clouds of gas, including gas 
in the vicinity of the progenitor, can leave absorption signa- 
tures on the spectra of SLSNe, and these may transmit the 
chemi cal composition of t hese otherwise un detectable sys- 
tems (jOuimbv et alJl201ll ; lBerger et al.ll2012h . If the SLSNe 
are, in fact, products of massive stars, then their volumetric 
rates will be entwined with th e cosmic star-forma tion history 
and the initial mass function |Tanaka et al.ll2012h . Finally, if 
the luminosities fall within a narrow enough range, or if they 
can be predicted from secondary indicators such as the light 
curve width, then SLSNe may serve as standard candles. 

Among the first clear SLSN discoveries are several con- 
tributions from the ROTSE-IIIb telescope at the McDonald 
Observatory. Although modest in size (the primary mirror 
is just 0.45 m), ROTSE-IIIb has a wide field of view (about 
1.85 x 1.85 degrees per exposure), which al lows large swaths 
of sky to be monitored l|Akerlof et al.f2 003L As the discover- 
ies are necessarily bright [M < 18.3), follow-up spectroscopy 
is relatively inexpensive. 

ROTSE- I IIb de tected one of the first SLSNe, SN 2006gy 
|Smith et all I2007I ). Classified as a Typelln from the 
narrow hydrogen emissio n lines apparent in its spectra 
llHarutvunvan et all 120061 : IPrieto et all 120061 : iFolev et all 
2006), SN 2006gy shone brighter than a typical Type la su- 
pernova at peak for 5 months and radiated over 10 51 erg 
in optical light alone. Other SLSN-II discove ries have fol- 
lowed including SN 2003m a (jRest et al.lboill ) SN 2008am 
dChatzopoulos et al.l l201lh . and SN2008fz (|Drake et all 
l2010h . Interaction with circum-stellar media clearly plays 
an important role in these events and may be responsible 
for their extreme luminosities (e.g. Smith fc McCravl 120071 : 



2011 



ISmith et al"1l2008bl : IChevalier fc Irwinl 

The nature of a prior ROTSE-IIIb discovery, 
SN2005ap, was not immediately apparent. The first spec- 
trum lacked the strong P-Cygni profiles typically seen in 
supernova spectra, but the rise in flux t oward shorter 
wave lengths did indicate at hot photosphere l|Quimbv et al.l 
120071 ). Later discoveries by the Paloma r Transient Factory 
" iRau et al.ll2009h and Pan-STARRS 
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l|Kaiser et alJl2010D of similar objects show that this class of 
objects is depleted in hydrogen and thus cannot produce its 
luminosity through interactions with a hydrogen rich CSM 
llPastorello et al. I l2010l ; lOuimbv et aTll201ll ; IChomiuk et al.l 
12Q11J). 

A rough estimation of the rate of ob jects similar to 
SN2005ap was presented in lQuimbvl (120081 ). The derivation 
used the discovery ratio of SN 2005ap to Type la supernovae 
in the ROTSE-IIIb sample (at that time) and a simplistic 
estimate of the search volume rati o for these to estim ate 
a rate relative to the Type la rateFl IMiller et al.l (I2009T ) fol- 



lowed this method with a larger sample from ROTSE-IIIb to 
derive an approximate rate for objects similar to SN 2005ap 
and SN 2008es relative to the Type la rate. They find about 
1 such event for every 350 Type la. They further note that 
the non-detection of su ch events from the KAIT SN search 
l|Filippenko et al.ll200ll ) implies that the rate in large galax- 
ies is less than 1/160 times the local Type II rate (note that 
2006gy was actually detected by the KAIT SN search in 
such a host, b ut it was not disco vered due to a selection 
bias) . Recently, ICooke et alj (|2012l ) estimated the high red- 
shift SLSN rate using two transients found at z — 2.05 and 
z = 3.90 in archival SNLS data. These were not spectroscop- 
ically confirmed, but they imply a SLSN rate at z = 2 — 4 
of ~ 400 Gpc^yr" 1 /^. 



lOuimbv et al.l (|2012h . have measured the volumetric 
rate of Type la supernovae using background discoveries 
from the ROTSE-IIIb searches. In this paper, we measure 
the rates of SLSNe from the same survey over the same pe- 
riod. It is not yet clear if SLSNe are produced from one 
or multiple channels, so we will supply the necessary pa- 
rameters for the reader to determine the approximate rates 
(or limits) of individual events from our survey, or to group 
events into physically related groups. We also supply rate 
measurements of the bulk SLSN-I and SLSN-II populations 
based on 1 and 3 discoveries, respectively. We further sup- 
ply the rate for an all inclusive category of SLSN-like events, 
which includes one additional (as yet unverified) SLSN can- 
didate in addition to the confirmed SLSN-I and SLSN-II. 
We discuss this event and the full sample in more detail in 

m 

These rates may be compared against prospective pro- 
genitor rates (including any high luminosity tail of nor- 
mal supernova) to help constrain the origin of these events. 
The rates of Type la supernovae and, in particular, core 
collapse events spanning the Type Iln through Type Ic 
classes, have previously been used to make su ch inferences 
and continue to be useful f o r such studies fseelGraur et all 
l201ll , iHoriuchi et"aH l201ll . iDahlen et ail I2OI2I . and refer- 
ences therein). 

For our rate estimation, we perform a Monte Carlo sim- 
ulation where the light curves of simulated SLSNe are com- 
pared against observations logged by ROTSE-IIIb to esti- 
mate the fraction of events we would recover over a variety 
of distances. We construct a sample of light curve templates 
from published SLSNe in H3.1l and use the peak magnitudes 
from these to estimate the pseudo-absolute magnitude dis- 
tribution (pAMD) , which is the intrinsic luminosity function 
mixed with the host absorption distribution, in >i3.2l The 
rates are presented in ^4j and discussion and conclusions 
are given in JSl 



2 THE ROTSE-IIIB SLSN SAMPLE 

Our sample is drawn from two similar surveys conducted 
with the RO TSE-IIIb teles cope: The Texas Supernova 
Search (TSS; lOuimbvl l2006h and the ROTSE Supernova 
Verification Project (RSVP; lYuanll2010P l. Both surveys used 
ROTSE-IIIb's wide field of view to canvas about 500 square 



We note that the volume ratio published in iQuimbvl {2008) 



is erroneous. Given the stated procedure, the calculation should 
have employed a volume ratio of about 27, not 3. 
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Table 1. Background SLSN-like Events Discovered with ROTSE-IIIb by Feb. 1, 2009. 



Name 


Disc. Date 


RA 


Dec 


z 


Type 


Reference 


2005ap 


Mar 3, 2005 


13:01:14.8 


+27:43:31 


0.283 


Ic 


Ouimbv et al. (2007) 


2006tf 


Dec 12, 2006 


12:46:15.8 


+11:25:56 


0.074 


iln 


Smith et al. f 2008a) 


2008am 


Jan 10, 2008 


12:28:36.3 


+15:34:49 


0.234 


Iln 


Chatzopoulos et al. (2011) 


2008es 


Apr 26, 2008 


11:56:49.1 


+54:27:26 


0.202 


11 


Gezari et al. (2009) 


Dougic 


Jan 21, 2009 


12:08:47.9 


+43:01:21 


0.191 


? 


Vinko et al. (in prep.) 



degrees on the sky with preference given to areas with high 
concentrations of galaxies in the local {D < 200 Mpc) uni- 
verse. 

Between November 1, 2004 and January 31, 2009, 
ROTSE-IIIb discovered 76 supernovae. All of these were 
spectroscopically confirmed by us or others in the commu- 
nity. The highest luminosity Type la fo und in our survey 
is SN 2007if at a bout M r = -20.4 mag ijScalzo et al.l 120101 : 
lYuanet al.lboioh . Here we consider only higher luminosity 
events, so we remove SN 2007if and 69 fainter objects. To cal- 
culate a volumetric rate that is unbiased by the local density 
enhancements (e.g. galaxy clusters) targeted by our search, 
we must further remove discoveries found about 50 Mpc be- 
hind the targeted objects or less. This cut eliminates what is 
perhaps th e best known SLSN discovery from ROTSE-IIIb: 
SN2006gy (|Smith et al.ll2007l : IOfek et al.ll2007f) . This object 
was found in a massive galaxy (NGC 1260) residing in the 
Perseus Galaxy Cluster, which we had specifically targeted. 

After these cuts, w e are left with the fi ve events listed 
in tableQ] SN 2005ap (|Quimbv et al.l 120071 ) was first de- 
tected before SN 2006gy, but it could not be identified as 
a SLSN until after the discovery of SN 2006gy reset the ac- 
cepted limits of the supernova peak absolute magnitude dis- 
tribution. SN2005ap is the first member of the hydr ogen- 
poor, SLSN -I group. SN2006tf (ISmith et~al] l2008al ) and 
SN 2008am (|Chatzopoulos et al.ll201ll ) both show hydrogen 
Balmer lines of relatively narrow widths in their spectra (in 
addition to broader components), and can be grouped in 
the SLSN-II categor y. The hydrogen features seen in the 
spectra of SN2008es (jMiller et al.ll2009l : iGezari et al.ll2009l ) 
lack these narrow emission peaks, and it was not until the 
initially hot, blue continuum cooled and faded after peak 
that broad P-Cygni profiles clearly emerged. This points to 
some differences in the pro genitor system-but not nec essar- 
ily profound differences (cf . iMoriva fc Tominagall2012T ) , and 
SN 2008es can be grouped with the other SLSN-II based 
on the eventual signs of hydrogen. Finally, ROTSE-IIIb 
detected a high luminosity event internally designated as 
"Dougie" (sometimes called ROTSE3 J120847.9+430121). 
Dougie exhibited a mostly featureless blue continuum, which 
was seen to redden over time; however, the broad features 
typical of supernovae were never observed. It is possible that 
Dougie is a tidal disruption event or an unusual AGN out- 
burst (although no x-ray emission was detected), but we do 
not rule out a connection to SLSNe at this time. Full details 
of this event will be presented later (Vinko et al. in prep.). 

In summary, we have just one SLSN-I (SN2005ap), 
three SLSN-II (SNe 2006tf, 2008es, and 2008am), and one 
additional, SLSN-like event (Dougie) to use in our rate cal- 
culations. If Dougie is a supernova, its lack of spectroscopic 



evidence for hydrogen would place it in the SLSN-I group, 
and we will supply SLSN-I rates with and without this event. 



3 LIGHT CURVES AND PSEUDO-ABSOLUTE 
MAGNITUDE DISTRIBUTIONS OF SLSNE 

A key factor in determining the rates of supernovae is ac- 
curately knowing the peak magnitude distribution with the 
effects of host absorption factored in. Here we devise repre- 
sentative light curves and pseudo-absolute magnitude distri- 
bution (pAMD) models in our unfiltered band pass for both 
SLSNe-I and SLSNe-II that can be used to determine our 
survey efficiency and thus rates. The pAMDs required for 
our rate study give the distribution of absolute magnitudes, 
without correction for absorption external to the Milky Way, 
that would be recovered by an ideal, volume limited survey. 
As the available SLSN sample has been gathered from non- 
ideal, flux limited surveys, we discuss how we can fit our 
volume limited pAMD models to the observed sample af- 
ter filtering the models for flux limited selection bias. There 
could always be additional bias, particularly the unknown 
publication bias, that we cannot account for. Given the small 
ROTSE-IIIb sample size, however, the precision of rates may 
yet be dominated by statistical errors, and if our pAMD 
models prove deficient in lower luminosity events, our rates 
will still stand as lower limits on the larger population. 



3.1 SLSN Light Curves 

We define a light curve sample using published photome- 
try of SL SNe. The sampl e consists mostly of the 18 events 
listed in iGal-Yaml |2012l ). The light curve for SN 1999as 
(|Knop et al.lll999l ) has not been published and we are thus 
unable to include it. The light curve sample includes the 
ROTSE-IIIb discoveries above. We have added unpublished 
observations of SN 20 06tf, Dougie, and a nother ROTSE-IIIb 
discovery, SN2010kd l| Vinko et alj|2010h . a SLSN-I that was 
discovered after the survey period considered here. Table[2] 
lists the SLSNe considered. 

We have gathered the observed photometry in each of 
the available pass bands, and we use these to estimate the 
rest frame absolute magnitudes in the ROTSE-IIIb band 
pass. When multiple bands are available, we select the sub- 
set of (sometimes several) pass bands that best overlap 
with the rest frame ROTSE-IIIb response. We convert the 
observed magnitudes to pseudo-absolute magnitudes with 
M = mx — Ax ~ kx, where m x is the observed mag- 
nitude in band X, n is the distance modulus, Ax is the 
Galactic extinction in band X, and kx is the k-correction 
term for converting the observer frame pass band into the 
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Table 2. Light curve sample of SLSN-like events. 



Name 



Type 



Abs. Mag™ Am4o 



Effective 
VoL-Time 6 



Effective 
Vol.-Time c 



Reference 



SN2003ma 


SLSN-II 


-21.58 


0.55 


0.0296 


0.0286 


SN2005ap 


SLSN-I 


-22.15 


1.36 d 


0.0503 


0.0288 


SCP06F6 


SLSN-I 


-22.11 


1.41 


0.0561 


0.0341 


SN2006gy 


SLSN-II 


-20.75 


0.64 


0.0113 


0.0349 


SN2006oz 


SLSN-I 


-21.67 d 








SN2006tf 


SLSN-II 


-20.53 


0.15 


0.0103 


0.0432 


SN2007bi 


SLSN-I 


-21.02 


0.35 


0.0182 


0.0450 


SN2008am 


SLSN-II 


-21.77 


0.24 


0.0653 


0.0503 


SN2008cs 


SLSN-II 


-22.02 


0.87 


0.0510 


0.0269 


SIN 2008fz 


CT CM TT 


—21.91 


0.70 


0.0494 


0.0311 


Dougie 


SLSN-I? 


-22.50 


2.66 


0.0332 


0.0129 


PTF09atu 


SLSN-I 


-21.59 


0.52 


0.0385 


0.0458 


SN2009jh 


SLSN-I 


-21.68 


0.69 


0.0397 


0.0416 


PTF09cnd 


SLSN-I 


-21.90 


0.55 


0.0565 


0.0452 


CSS100217 


SLSN-II? 


-22.79 


0.18 


0.2870 


0.0546 


SN2010gx 


SLSN-I 


-21.54 


1.54 


0.0220 


0.0288 


PSf-lOky 


SLSN-I 


-21.92 


1.15 


0.0433 


0.0317 


PSl-lOawh 


SLSN-I 


-21.90 




0.0356 


0.0273 


SN2010kd 


SLSN-I 


-21.08 


0.28 


0.0245 


0.0562 



Rest et al. (20111 
Ouimbv et al. (2007) 
Barbarv et al. (20091 
Smith et al. (2007) 
Leloudas et al. (2012) 
Smith et al. (2008a) 
Gal- Yam et al. (2009) 
Chatzopoutos_et_aL (2011) 
Gezari et al. (2009) : Miller et al. (2009) 
Drake et al. (2010) 
Vinko et al. (in prep.) 
Quimby et al. (2011) 
Quimby et al. (2011) 
Quimby et al. (2011) 
Drake et al. (2011) 
Pastorello et al. (2010) 
2011) 
2011) 



Chomiuketal 
Chomiuk et al 
Vinko et al. (in prep.) 
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"Pseudo-absolute magnitudes are in the unfiltered ROTSE-IIIb rest frame system and not corrected for host absorption. 
1 Effective volume-time in units of Gpc 3 yr KZ^ using only the object's light curve and an assumed intrinsic Gaussian distribution of 

peak magnitudes (cr = 0.3 mag ). 

Effective volume-time in units of Gpc 3 yrhZ-j* using only the object's light curve scaled to the average peak magnitude of the group 

(-21.7 ±0.4 mag for SLSN-I and -21.4 ±0.6 mag for the SLSN-II). 
^Estimate based on extrapolation of the light curve. 



rest frame ROTSE-IIIb band. We do not include a term to 
correct for host absorption, hence our values are "pseudo- 
absolute". Our unfiltered ROTSE-IIIb magnitudes are cali- 
brated against the USNO-B1.0 R2 magnitudes, and we in- 
clude the conversion from observed magnitudes in the AB 
system to the Vega system in the k-correction term when 
applicable. 

To properly estimate the k-corrections, we need to know 
the intrinsic spectral energy distribution (SED) of each tar- 
get as a function of time. As this information is not always 
available, we make the following approximations. 

For the SLSNe-I, we assume an augmented, cool- 
ing black body for the SED. Using the analysis of 
IChomiuk et"ail(|201ll ). we fit a cooling curve to the sample of 
SLSN-I with spectroscopically derived black body temper- 
atures. We add the best fit black body temperature mea- 
sured from the spectra of SN 2007bi about 48 rest frame 
days after peak, which fits in with the overall cooling trend. 
To approximate the effects of metal absorpt ion in the UV 
l|Pastorello et alj|201ol ; IChomiuk et al.ll201ll ). we scale the 
template flux below 3000 A by a linear function that grows 
from at 912 A to 1 at 3000 A. We compare the observed 
spectra of several SLSN-I covering a range of rest frame 
wavelengths to our assumed templates at similar phases in 
Figure[T] As we are ultimately interested in the response over 
ROTSE-IIIb's very broad, unfiltered band pass (or the com- 
bination of multiple broad bands) , we assume that the effects 
of individual line features missing from the templates to be 
minor. For example, if we compare to the k-corrections de- 
rived from the actual spectrum of PTF09cnd shown in Fig- 
ure [T] to those derived with our SED model, then, assuming 
a redshift of z = 0.17, our SED model (at the appropriate 



100.00 




2000 4000 8000 

Rest Wavelength (A) 

Figure 1. SLSN-I spectral templates (thick blue, green, and 
red lines) compared to the actual spectra of supernovae at 3 dif- 
ferent phases (20 days before maximum, near maximum, and 27 
days after maximum in the rest frame ROTSE-IIIb system). The 
spectra of PTF09cn d, SCP06F6, and SN2 010g x are taken from 
lOuimbv et~aH [|201ll) , [Barbarv et al.l (2009), and iPastorello et alj 
l|201fj ). respectively. The yellow shaded area marks the ROTSE- 
IIIb rest frame band pass. 



phase) gives k-corrections in the g band that are ~ 0.1 mag 
too faint, accurate to less than 0.01 mag in the r band, and 
less than 0.1 mag too bright in the i band. Averaging the 
three bands should thus result in reasonably accurate abso- 
lute magnitudes (we will adopt a 0.2 mag offset to estimate 
systematic error in our rates). 

Near maximum light, which is the most important phase 
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Figure 2. Blackbody temperature of SLSNe-II as a function of 
the parameter, rj, described in the text. T he symbols with er - 
ror bars plot indivi d ual measurem e nts frq mlDrake et alj J2 01ll). 
ISmith et al l 1120071), ISmith et alj ll2008ah . IChatzopouIos et alj 
(|201lh . and lMiller et al.1 l|2009l ). The black line is the best fit cubic 
polynomial for —0.4 < r\ < 0.8 and a constant above and below 
these limits. 



for our rate calculation, the features present in the rest frame 
ROTSE-IIIb are very weak in comparison to normal super- 
novae, so the temperat ure of the black body is the main 
concern. The results of IChomiuk et all |201ll ) suggest that 
SLSN-I follow similar cooling curves, although there may be 
a range of peak temperatures. We note that the SDSS-II 
photometry of SN 2006oz indicat es a nearly constant p hoto- 
spheric temperature with time (jLeloudas et al.l 12012). and 
the earliest phases are cooler than our simple model. The 
peak of this particular event is not well constrained, so we 
cannot include it in our efficiency calculations. 

For the SLSNe-II, we assume a simple Planck function 
for the SED. This is motivated by events like SN2008es, 
which was well observed and found to be we ll fit by a 
black body (|Miller et al.l2009l;lGezari et al. 120091 ). and other 
SLSNe-II like SN 2006tf (ISmith et al.l2008al ) and SN 2008am 
l|Chatzopoulos et al.l 1201 ll ). which, apart from their hy- 
drogen emission lines, are reasonably well approximated 
by black bodies. The temperature measured at peak for 
SN2008es (~14,000K) is much higher than for SN2006tf 
(~8000 K) , but both objects cool in a similar manner as they 
fade. We find that the temperature evolution of SLSNe-II 
can be parameterized as T(r\) — 11.8 — 15.37) + 14. 4^ 2 —4.2rj 3 
(in units of 10 3 K; see Fig. [2}. The parameter, 77, is a 
function of the peak absolute magnitude in the rest frame 
ROTSE-IIIb band, M peak , and the decline from peak, AM, 
such that 77 = 0.25(Af pe afe + 22) — AM before peak and 
77 = 0.25(M P eafe + 22) + AM after. The 77 parameter thus 
corresponds to the drop in magnitudes below peak (with 
negative values for pre-maximum phases) normalized to the 
behavior of a M = —22 source. We limit 77 to —0.4 < 77 < 0.8, 
so the temperatures can be up to 20000 K early on, but never 
below 6700 K at later times. We have also experimented with 
other ad hoc parameterizations of the temperature evolu- 
tion, but these made no significant changes to our results. 

To construct the final light curve templates, we fit low 
order polynomials to the photometry (after conversion to 



the rest frame ROTSE-IIIb band pass as described above). 
In some cases, we smoothly combine fits to the rising and 
declining phases to better capture asymmetries in the light 
curves without resorting to higher order polynomials. Since 
the date of maximum light depends on the k-correction, and 
the k-corrections depend on the phase relative to maximum 
light (and also the peak magnitude in the case of SLSNe-II), 
we begin with an initial guess for the date of maximum and 
iterate the k-corrections and polynomial fits until they con- 
verge. The resulting templates are shown in Figure [3] Note 
that our light curve templates are constrained by observa- 
tions taken, on average, every 3.5 rest frame days within 20 
days of maximum light (without double counting the multi- 
ple bands usually obtained on the same night). 

Several of the supernovae in our light curve sample 
have both unfiltered ROTSE-IIIb photometry and filtered 
observations. After k-corrections the ROTSE-IIIb measure- 
ments are about 0.1 mag fainter than the (similarly cor- 
rected) filtered photometry of SN 2006gy and SN 2008es, and 
the ROTSE-IIIb measurements are about 0.1 mag brighter 
for Dougie and SN 2010kd. The largest discrepancy seen is 
for SN 2006tf, for which the ROTSE-IIIb measurements are 
about 0.2 5 mag brighter than the filtered photometry avail- 
able from ISmith et al.1 (|2008a) after k-corrections. The ob- 
served ROTSE-IIIb magnitudes agree with the observed R- 
band measurements to within 0.1 mag over the same range 
(out to about 100 days after maximum light), so the dis- 
crepancy in k-corrected, absolute magnitudes is probably 
due in part to the k-corrections themselves (as noted above, 
we will assume a 0.2 mag systematic error in magnitudes for 
the rate measurements). In the case of SN2006tf, we choose 
to shift our unfiltered photometry to match the filtered mea- 
surements. This preserves the constraints on maximum light 
from the ROTSE-IIIb observations, which precede the fil- 
tered data by up to 1 month. 

CSS100217 ijDrake et al.1 l201ll ) appears to be dis- 
tinctly brighter than the rest of the objects considered, 
and, integrating its light curve, it radiated far more en- 
ergy than the others over the period shown in Figure [3] 
It was also located precisely coinci dent with a know n 
AGN_(SDSS J102912. 58+404219. 7; cf. I Greene fe~Holl2007l ). 
iDrake et al.l ()201ll ) favor classifying this object as a SLSN- 
II because its variability is larger than normal AGN flares, 
but given that its amplitude is also larger than even the ex- 
treme supernovae considered here, a similar argument could 
be used to disfavor classification as a supernova (see fur- 
ther discussion in SJ5J. Similarly, Dougie is brighter than the 
other remaining events, it rises and fades faster, and, as 
noted above, it shows a spectral evolution that is distinct 
from most SLSNe. While each of these may yet prove to be 
supernovae, the evidence against this interpretation is con- 
siderable. We will thus calculate the SLSN rates below with 
and without these questionable events. 



3.2 SLSN Pseudo-Absolute Magnitude 
Distribution 

From the light curve templates constructed above, we derive 
the approximate pseudo-absolute magnitude distribution of 
SLSNe-I and SLSNe-II. Ideally, one would like to have a vol- 
ume limited sample from which to derive these distributions. 
This is simply not available, and we are further limited by 
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Figure 3. Light curves of SLSN-like events in the rest frame 
unfiltered ROTSE-IIIb band pass. Blue curves are for SLSN-I, 
and red curves mark SLSN-II. The green curves are for two objects 
of uncertain origin. 



the small numbers in the flux limited samples. Given these 
limitations, it may be tempting to simply take the peak 
magnitudes from the light curve sample, measure the av- 
erage and standard deviation, and then assume a Gaussian 
distribution with these parameters (and we do consider such 
distributions for the rate calculations below for reference). 
However, this will generally be problematic since the ob- 
served sample is flux limited and thus lacking representa- 
tion from the lower luminosity members of the true popu- 
lation. The mean and standard deviation from the observed 
samples will thus be biased. Additionally, real, astrophysical 
sources that lie in star-forming galaxies will be obscured by 
dust. This will skew an intrinsically Gaussian distribution 
and produce an extended faint tail to the pAMD. 

We thus attempt to derive realistic pAMDs from the 
data available by assuming the intrinsic populations obey a 
Gaussian distribution for simplicity and including the skew- 
ing effects of dust. To translate these ideal models into the 
distributions recovered by a flux limited search, we simply 
weight the ideal distributions by the relative volumes ex- 
pected from such a survey. We can then find the best pAMD 
(the one that would be found from an ideal, volume limited 
survey) by varying the input model until the best match is 
found between the weighted model and the observed sample. 
As the model separates the effects of dust from the intrinsic 
scatter, insights into the later may be gained if the correct 
dust model is used. SN 2006oz was not observed at maxi- 
mum light, so we must exclude it from what follows. 

The sample is based on discoveries drawn from several 
different surveys, each with their own selection biases. The 
ultimate selection function for our light curve sample is thus 
complex and largely unknown. To simplify matters, we as- 
sume the aggregate population is selected from a hypothet- 
ical, flux limited survey. This should be a reasonable as- 



sumption since the contributing surveys are flux limited (not 
targeting specific luminosity ranges or host galaxy types) 
at their core, with spectroscopic follow-up adding a second 
roughly flux limited layer to the discovery process. There 
could be additional biases, however, such as a preference for 
publication of the higher luminosity events, that may further 
skew our light curve sample. Such effects should typically 
limit the inclusion of lower-lumi nosity events (i n particular, 
we note that the sample from iGal-Yaml |2012| is based on 
events that peak above —21 mag in any given photometric 
system). This could lead us to overestimate our survey effi- 
ciency and thus underestimate the rates. 

We consider simple models for the pAMDs based on 
intrinsic luminosity functions that follow a Gaussian distri- 
bution combined with host absorption drawn from an ex- 



ponential distribution of the form P{Av) oc e 



-A v Jr 



(the 



V-band is the best match to the ROTSE-IIIb band when 
redshifted to z ~ 0.2). We adopt r = 0.6 for the SLSN- 
I and SLSN-II populations because these events are very 
likely to be connected with massive stars and the absorp- 
tion distribution expected for objects originating from such 
stars is approximately fit with this exponential distribu- 
tion (|Hatano. Branch fc Deatonlll99cf ). As shown below, this 
choice also provides a reasonable match to the observed dis- 
tributions of SLSN in the light curve sample after transform- 
ing the pAMD models to the distributions expected from a 
flux limited survey. 

With r fixed, we vary the model's peak (intrinsic) mag- 
nitude, M, and its Gaussian width, a to find the best match 
to the light curve sample. We measure the maximum offsets 
in the cumulative distributions of the models and data and 
we record the Kolmogorov-Smirnov probability that the ob- 
served distribution was drawn from each model. Given the 
limited SLSN-I and SLSN-II sample sizes, a variety of values 
are plausible, but we find the best agreement for the SLSN-I 
sample around M = -22.0 and a = 0.3, while the SLSN-II 
population is best matched with values around M = —21.4 
and a = 0.6. The peaks of the intrinsic distributions can- 
not be much brighter (e.g. M > -22.3 for the SLSN-I), but 
fainter values are plausible if the widths of the distributions 
are larger. We did not allow the absorption distribution to 
vary, but smaller values of r would decrease the expected 
number of the lower luminosity events, which is less pre- 
ferred (but not excluded) . We considered observed distribu- 
tions both with and without CSS100217 and Dougie, but 
the differences are negligible. The resulting pAMD models 
are shown in Figure [4] 

We have included SN2007W and SN2010kd in the 
SLSN-I light curve sample even though their light curves 
decline (and likely rise) much more slowly than the other 
SLSN-I, which could point to a different physical origin 
|Gal-Yam et al.ll2009l have concluded that SN 2007bi is a ra- 
dioactively powered, pair-creation supernova, but this pos- 
sibility is excluded for others in the SLSN-I sample). Even 
with these events, the width of the intrinsic luminosity func- 
tion could be rather narrow. If there is a publication bias 
against lower luminosity SLSN-I, then the true distribution 
could be broader than the available evidence suggests. It is 
also possible to explain the observed scatter in SLSN-I peak 
magnitudes from the effects of dust alone, with zero intrinsic 
dispersion (assuming the dust model adopted above). 

We estimate the uncertainty in the peaks of our distri- 
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Figure 4. Peak pseudo-absolute magnitudes for SLSN-I (blue) 
and SLSN-II (dashed red) in the light curve sample. SN 2006oz is 
not included because its peak magnitude is not well constrained. 
The locations of CSS100217 and Dougie, which may or may not 
actually be SLSNe, are marked with the green arrows on the top 
plot (cumulative distribution) and dashed-dotted green boxes on 
the lower plot (binned distribution). The dark blue and dotted 
pink curves show the models used in the survey efficiency calcu- 
lations adjusted for the selection effects of a flux limited survey 
for comparison to the observed sample. The SLSN-I model has an 
intrinsic Gaussian width of just a = 0.3 mag. The skewed faint 
end tail of the model is a result of the assumed host absorption 
distribution. 

buttons by varying the temperatures assumed in calculating 
the k-corrections for our light curve sample. To determine 
the systematic error, we adjust the model temperatures by 
±1000 K, and add to this an estimate of the statistical un- 
certainty in the peak magnitudes. 



4 THE SLSN RATES 

In this section, we combine the ROTSE-IIIb sample pre- 
sented in section [5] and the light curves and pAMDs infered 
from the published sample of SLSNe (presented in section 
[5J to derive the volumetric rates of SLSN-I, SLSN-II, and 
all SLSN-like objects (which may include additional, rare 
phenomena that may or may not be true supernovae). The 
rates are calculated from, 



with N the number of events detected, V5 the co-moving 
volume element for the i th distance bin, e; the corresponding 
survey efficiency for that bin, and ti the proper time for the 
survey in each bin. We refer to the denominator in eqn. 
[T] as the "effective volume-time" in what follows. We first 
describe the Monte Carlo simulations employed to calculate 
the survey efficiency before presenting the actual rates. 



4.1 Monte Carlo Simulations 

Following the procedure described in IQuimbv et al.1 l|2012T ). 
we have run Monte Carlo simulations to determine how ef- 
ficient our ROTSE-IIIb search has been in selecting SLSNe. 
We define our survey efficiency, a, as the fraction of SLSNe 
in a given volume that we are likely to select. There are 
two main factors in selecting SLSNe: 1) are the observations 
deep enough to detect a given SLSNe at a given distance, 
and 2) did we observe the appropriate sky location during 
the right time to catch the transient during outburst (rel- 
evant only if the observations are deep enough). The first 
factor results from the flux limited (i.e. not volume limited) 
nature of our survey. At a given distance we may detect some 
of the more luminous SLSNe but miss some of the fainter 
events. The second factor is a result of our sampling func- 
tion (or "cadence"). Due to weather, observing seasons, and 
time constraints, we have not observed the same fields every 
day, and some fields have far less covera ge that others (see 
figures 1 and 2 from IQuimbv et ai1l2012r ). so we could have 
missed even bright SLSNe in coverage gaps. We thus per- 
form a Monte Carlo simulation in which we compare SLSNe 
light curves at various sky locations, distances, and random 
explosion dates to the actual survey data to determine what 
fraction of simulated events we are likely to have selected. 

We determine our selection efficiency in distance bins 
from 40 to 4000 Mpc (beyond this our selection prob- 
ability is zero, as we will show below). For each dis- 
tance bin, we convert the absolute magnitude light curves 
from section [3] to their expected (or "observed") magni- 
tudes by adding the distance modulus, the same crude k- 
corr ections described above, and the G alactic extinction 
from lSchlegel. Finkbeiner fc D avis ( 1998) for a given sky po- 
sition. Next we randomly select a date for maximum light 
between Nov 1, 2004 and January 31, 2009. The observed 
magnitude can then be computed for any day over the survey 
period. We compare the expected magnitudes for our sim- 
ulated SLSNe to the actual data recorded to determine the 
fraction of events likely to have been selecte d. To do this, 
we us e the detection efficiency curves from IQuimbv et al.l 
(2012), which give the probability of detecting an object as 
a function of magnitude. These detection efficiency curves 
take into account the seeing quality for the image as well as 
the "limiting magnitude." 

In the actual survey, the selection of targets was done by 
humans reviewing an automatically generated "short list" of 
candidates that was constructed from targets with at least 
a 5cr detection in a nightly image stack, and at least 2.5a 
detections in the first and second halves of the night (in 
addition to ot her cuts on shape pa rameters and increase in 
brightness; see IQuimbv et all [20 12 for full details). Thus, in 
our simulations we calculate the detection probabilities on 
each of the three stacks, draw three random numbers (uni- 
formly distributed between and 1), and we only count a 
simulated source as recovered if each of the three probabili- 
ties is larger than their respective random draws. 

As in IQuimbv et all l|2012t) , we compare the expected 
distributions in the distance, observed peak magnitudes, 
pseudo absolute magnitudes, and the number of nights a 
simulated source is detected from the simulations to the ob- 
served data. Given our small sample size, we lack the statis- 
tical power to identify differences between how the simulated 
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sources are selected versus the actual survey data; however, 
we do note that not one event in our observed sample was 
detected on fewer than 5 nights. In contrast, the simulations 
typically predict that ~60% of o ur discoveries should only be 
detectable on 4 or fewer nights. iQuimbv et all (|2012h found 
a similar disagreement between the models and data for a 
larger sample of Type la supernovae, which could be resolved 
by considering only the events (real or simulated) detected 
on 5 or more nights. The dearth of objects detected on just 
1-4 nights is presumably attributed to human bias in the se- 
lection of targets. We assume a similar effect for our SLSN 
sample and limit the models and data to events detected on 
5 or more nights. With this cut, the four model distributions 
considered appear to be well matched with the observations. 
If we were to ignore this selection cut, the rates found below 
would be lowered by a factor of 2.5. 

4.2 SLSN-I Rates 

We have only one discovery, SN 2005ap, from which to de- 
rive the SLSN-I rate. Obviously, our rate will be subject to 
the eccentric properties of small number statistics, and care 
should be taken in interpreting these results. 

We first perform Monte Carlo simulations using the 
SLSN-I light curves from section [3] (excluding Dougie for 
the moment) and the the SLSN-I pAMD model derived in 
section 13.21 to determine our effective volume-time. Under 
these assumptions, the effective volume-time of the ROTSE- 
Illb search is 0.0315 Gpc 3 yr 1 fe^ 3 , which gives a rate of 
(32^2g) events Gpc -3 yr -1 h?i, where the err or bars accoun t 
for only the statistical Poisson fluctuations (Gehrels 1986). 

To estimate our systematic error, we recalculate our 
survey efficiencies with the input pAMDs shifted 0.2 mag 
brighter or fainter. This is to account for possible systematic 
error in our k-corrections, which may bias the peak magni- 
tude distributions. Given our small sample size, the statis- 
tical errors dominate the systematics, which we find to be 
+7° events Gpc" 3 yr" 1 hf-i f° r tne SLSN-I. We note, however, 
that due to selection (and particularly publication) bias, our 
pAMD may not accurately represent the full population of 
events physically connected to the SLSNe-I. In this case, our 
rates are valid only for the luminosity range considered and 
are simply lower limits for the greater population. 

Figure [5] shows our survey efficiency as a function of 
luminosity distance. The expected distance distribution of 
SLSN-I is also shown with the solid blue curve. This latter 
curve is calculated by multiplying the efficiency curve by 
the relative volume and proper time in each distance bin. 
Although the sensitivity curve drops off above ~300 Mpc, 
this is initially compensated for by the nearly cubic rise in 
volume with distance. Due to this rise in volume, we are 
sensitive to SLSN-I out to more than 2000 Mpc even though 
we may detect less than 1% of this population. We expect 
to find about 97% of our SLSN-I between 200 and 2000 Mpc 
(roughly 0.05 < z < 0.4), and half should be farther than 
about 800 Mpc. Note that at D L ~ 1440 Mpc, SN2005ap 
falls comfortably with in our expectations. The average ef- 
fective redshift for our SLSN-I rate measurement is thus 
2 = 0.17. 

As a check, we also provide effective volume-times as- 
suming only the light curves of individual supernovae and 
two separate assumptions on the peak magnitude distribu- 



tee — i — ' — ' ' ' i 1 ' — ' — i — i — ' ' ' i <- 

'- SLSN-I pAMD 

: SLSN-II pAMD 




Luminosity Distance (Mpc) 

Figure 5. Survey efficiencies (dashed lines) and the relative 
distance distributions (solid lines) of SLSN-I (blue curves) and 
SLSN-II (red curves) expected from these. The luminosity dis- 
tances for the SLSNe in the ROTSE-IIIb sample are indicated 
with arrows. 

tion in table [2] The values in the fifth column assume that 
the distribution is a simple Gaussian (no skewing by host ab- 
sorption) that takes the reported peak magnitude of the su- 
pernova for the peak of the distribution and assumes a width 
of a = 0.3 mag. This width was chosen to represent the (as- 
sumed) intrinsic scatter in events otherwise identical to each 
object listed and is based on the intrinsic scatter of SLSN-I 
estimated above. For SN2005ap, the effective volume-time 
is 0.0504 Gpc 3 yr 1 ft^ 3 under these assumptions. From equa- 
tion [T] this leads to a rate of (20lt*|) events Gpc" 3 yr" 1 
(statistical error only). 

The sixth column in table [2] gives the effective volume- 
times again assuming only the single light curve of the su- 
pernova and a simple, uncorrected Gaussian distribution 
but with the peak and width set by the average and stan- 
dard deviation of the apparent SLSN-I sample, respectively 
(M = —21.7 and a = 0.4 mag). With these assumptions, the 
effective volume-time drops to 0.0288 Gpc 3 yr 1 h^i- The rate 
for this distribution is then (351*9 ) events Gpc" 3 yr" 1 h^. 

Our original measurement is fully consistent with these 
two later assessments. Again, the small sample size (one ob- 
ject) results in large statistical errors which dominate sys- 
tematic concerns. As the procedure used to estimate the 
pAMD model is generally less prone to bias than the later, 
simple Gaussian assumptions, we retain its effective volume- 
time for our best measurement of the SLSN-I rate. 

If we count Dougie as a SLSN-I, then the effec- 
tive volume-time (including Dougie and all SLSN-I light 
curves in the Monte Carlo simulations) decreases slightly to 
0.0291 Gpc 3 yr 1 . The combined rate for these two events 
is then (68I44) events Gpc" 3 yr" 1 ft 3 ! (statistical error only). 

4.3 SLSN-II Rates 

Our ROTSE-IIIb sample includes three events with clear sig- 
natures of hydrogen in their spectra: SNe 2006tf, 2008am, 
and 2008es. Although they share some similar properties, 
there are differences as well. In particular, the available ob- 
servations of SN 2008es do not show the obvious narrow 
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emission lines from hydrogen. These are present in the spec- 
tra of SNe 2006tf and 2008am and suggest an interaction of 
the supernova ejecta with slow moving circum-stellar ma- 
terial as a contributing power source. It could be, however, 
that SN 2008es experience d a similar interaction prio r to the 
first available spectra (cf. iMoriva fc Tominagall2012l ). 

Given the apparent differences, one could follow a simi- 
lar procedure as outlined above and calculate three separate 
rates using the effective volume-times in Table [2] Here we 
use the full ensemble of SLSN-II light curves available in the 
light curve sample and the pAMD model from section [3721 to 
place limits on the bulk, hydrogen rich SLSN population. 

Using the SLSN-II pAMD model, we find an ef- 
fective volume-time of 0.0199 Gpc 3 yr 1 (this would 
be 6% larger if the light curve of CSSI002I7 were in- 
cluded). From equation [1] this gives a SLSN-II rate 
of (ISI^gj 1 ) events Gpc -3 yr -1 (statistical error only) 
for our sample of three events. Performing the same 
check described above, we estimate a systematic error of 
^33 events Gpc -3 yr -1 . Again, we stress that there could 
be lower luminosity events that are physically related to 
those considered in determining the pAMD. In this case, 
our rate is a lower limit on the greater population. 

From Figure [5] it is apparent that like the SLSN-I, we 
are sensitive to SLSN-II in roughly the 200 to 2000 Mpc 
range. The average redshift for our SLSN-II rate measure- 
ment is z = 0.15. 



4.4 SLSN-Like Rates 

Finally, we measure the total rate of SLSN-like events in- 
cluding possible tidal disruption events or rare, supernova- 
like outbursts from (uncatalogued) AGN. Our survey was 
not intentionally biased against events located near galaxy 
nuclei (for example, SN 2006gy was found less than 0.3 pixels 
from the core if its host), but we did exclude a few sources 
detected coincident with known AGN. Our detection effi- 
ciency was also calculated for blank sky locations only, so it 
may be an overestimation for sources on top of bright galac- 
tic nuclei. So once again, our rate may be considered a lower 
limit. 

We include Dougie with the other 4 background events 
and assume the effective volume-time is simply given by 
the average of the SLSN-I and SLSN-II values found above 
(including Dougie and CSSI002I7 in the Monte Carlo 
simulations). Since the pAMDs for these two groups are 
broadly similar, this should be a reasonable approxima- 
tion. With these assumptions, the total SLSN-like rate is 
(199*^1^) events Gpc -3 yr -1 hfi, where the first and sec- 
ond errors quoted are statistical and systematic, respec- 
tively. 



5 CONCLUSIONS AND DISCUSSION 

We have produced a set of SLSN light curve templates cor- 
rected to the rest frame, unfiltered ROTSE-IIIb band pass 
and used these to estimate the pseudo-absolute magnitude 
distributions of SLSN-I and SLSN-II. With these and our 
sample of ROTSE-IIIb discovered SLSN (and a SLSN-like 
event), we have calculated the volumetric rates in the local 
volume (z ~ 0.2). Compared to the z ~ 0.2 core collapse 



supernova (CCSN) rate (cf. iBotticella et al.llioosh '). we find 
that there is about one SLSN-I or SLSN-II event for every 
400 to 1300 CCSNe, or one SLSN-I for every 1000 to 2 x I0 4 
CCSNe (these fractions may be halved if t he full popula- 
tion of faint or obscured CCSNe is included; iHoriuchi et al.l 
l201ll ). The total SLSN-like rate is also similar to the lo- 
cal rate of sub-energetic, long -duration gamma-ray bursts 
(LGRB; ISoderberg et al.ll2006h . although the SLSN-like and 
LGRB samples are currently too small for a definitive com- 
parison. 

If the SLSNe rates simply track the cosmic star for- 
mation rate ( CSFR ), then, using the parameterization of 
lYiiksel et ail (|2008T l. the CSFR and thus SLSN rate should 
be about 5 times larger at z ~ 3 than we have measured for 
the nearby ROTSE-IIIb sample. Given the large errors, this 
could b e roughly consisten t with the rate of SLSNe at z ~ 2 
an d 4 llCooke et al.ll2012l L However, the search performed 
bv lCooke et all (|2012l ~ was likely biased to longer duration 
events, and they suggest their discoveries represent only the 
radioactively powered subset of SLSNe (and only a fraction 
of those). In this case, their high redshift discoveries could 
imply an enhancement over a simple scaling by the CSFR 
when compared to an appropriate fraction of our low red- 
shift rates, as they have noted. 

The SLSN-I in the light curve sample are particularly 
interesting. Although these events were observed in a variety 
of (rest frame) band passes and we apply only a simplistic 
k-correction, they still appear to constitute a surprisingly 
uniform class when transformed to the ROTSE-IIIb system 
(see Fig. EJ) . The light curves show some quantitative differ- 
ences but tend to rise smoothly over around 50 rest frame 
days, which is considerably longer than typical Type la su- 
pernovae. Several SLSN-I also show smooth declines that 
almost mirror their pre-maximum behavior. The peak mag- 
nitudes of the SLSN-I in the light curve sample have an aver- 
age of —21.7 mag and a standard deviation of ±0.4 mag. Our 
best estimate of the intrinsic scatter (after accounting for 
absorption in the host galaxies and selection bias) is about 
±0.3 mag, but this is not well constrained. These values are 
based on a small sample (10 events) and change negligibly 
if Dougie is added. Given that the peak magnitudes of the 
SLSN-I sample are relatively tightly distributed, it is worth 
considering if this is a mere selection bias or indicative of a 
population that may prove useful as standard candles. 

First, we consider the high luminosity cutoff to the dis- 
tribution. The published sample of SLSN-I shows a cap 
in peak absolute magnitudes close to M ~ —22. If there 
are such events with significantly brighter peaks they must 
be exceedingly rare or they would dominate flux limited 
searches. It is thus reasonable to conclude that most SLSN-I 
will peak around —22 mag or fainter. 

Next, we consider the lower luminosity limit for the 
SLSN-I available in the light curve sample. Formally, we 
have used an arbitrary cut on peak luminosities (M poa k < 
—21) to select the SLSN-I for our light curve sample, but we 
note that we could use other, luminosity independent criteria 
to define this sample. In particular, the SLSN-I in the light 
curve sample contain the only supernovae we know of that 
show Oil features in their maximum light spectra. Other 
supernovae with more normal peak luminosities may also 
display such features at very early times, but these quickly 
disappear as the ejecta cool and oxygen recombines (cf. 
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the s pectral series of SN 2008D presented by iModiaz et alj 
2009). As a practical matter, this definition is perhaps less 
useful than a simple luminosity cut as it requires reasonably 
high signal to noise spectra taken at the right phase and 
with the proper wavelength coverage to catch the relatively 
weak Oil features in the 4400 to 5800 A range. Indeed, The 
wavelength coverage and signal to noise of SCP06F6 and 
PSl-lOky preclude a definitive assessment of such features 
at maximum light, but based on the available spectral infor- 
mation (particularly the high temperatures implied and the 
similar features at 2100 to 2700 A), we find this to be likely. 
There are no spectra of SN 2007bi near maximum light, but 
we suspect the Oil features would have been visible at max- 
imum as well considering later phase spectroscopic matches 
to the SLSN-I sample. 

There are no known supernovae of normal luminosities 
that show Oil at maximum light (events with this feature at 
maximum will be called SN 2005ap-like hereafter) , and this 
signature may uniquely identify a potentially useful pop- 
ulation of high luminosity supernovae. Based on the con- 
siderations below, we find it likely that all SN 2005ap-like 
events will be superluminous because the Oil lines require 
high temperatures and the long expansion of the ejecta up 
to maximum light necessitate a large radius. The maximum 
light spectra of SLSN-I are reasonably well approximated by 
a black body (at leas t in the optical range where line blan- 
ketin g is not an issue; iPastorello et all 2010, C homiuk et al.l 
l201ll ; see Fig. [[]), so the peak luminosity of SLSNe-I will be 
approximately, L = 4irR 2 aT 4 (with a small dilution term to 
correc t for the depa r ture from a black body, which we ne- 
glect ) . |Hatano^eFj2 (| 19991 ) have shown that Oil features are 
temperature sensitive. These lines should be replaced by the 
usual OI lines for temperatures below about 12000 K, and by 
OIII lines for temperatures above 15000 K. The typical ex- 
pansion velocity of supernovae, 10 4 kms -1 over a one month 
period gives a radius of ~ 3 x 10 15 cm. Thus any supernova 
that remains hot enough to display prominent Oil features 
after a one month rise to maximum light should have a peak 
absolute magnitude brighter than about M — —21. 

More detailed theoretical modeling, beyond the scope 
of this paper, is required to verify this simplified picture 
and determine if all SN 2005ap-like supernovae must nec- 
essarily reach such high luminosities. If the observed floor 
to the peak magnitude distribution proves to be physical 
and not just a selection effect, then the peak magnitudes of 
SN 2005ap-like supernovae are confined to a relatively nar- 
row range. This gives potential for SN 2005ap-like super- 
novae to one day serve as standard candles visible as far as 
L* galaxies can be seen and clearly warrants further study. 

We next consider if SLSN-I obey a luminosity-decline 
relation analogous to the Phillips relation observed for the 
standard bearer of standard candles, Type la supernovae. 
In Figure [6] we plot the change in magnitude between peak 
and 40 days after peak, A77140, against peak pseudo-absolute 
magnitude. SLSNe may be better differentiated in their de- 
clines from peak at even later phases, but the observations 
are harder to come by. Errors shown are the estimated sta- 
tistical error plus an estimate of the systematic error based 
on changing the temperature of the SED used in the k- 
correction, as discussed above. 

Using the Baye s ian linear regression routine, 
linmix_err .pro (|Kellvl [2007), we do not find a signifi- 
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Figure 6. Photometric decline in magnitudes between peak and 
40 days after peak for the SLSN light curve sample. The open 
orange diamond shows where SN 2006gy would lie if corrected 
for host absorption. The open circle marks SN 2005ap, which was 
not observed on day +40 or after; the A7TI40 shown for SN 2005ap 
is an extrapolation of the earlier light curve. Errors include the 
estimated systematics. 



cant correlation for 9 SLSN-I plotted in Figure [6] If we 
include Dougie, however, there may be a weak correlation 
between the peak magnitudes and decline parameter in the 
sense that brighter events have faster declines. This may 
simply be a selection effect, however, which would tend to 
remove the fainter, faster declining events. 

Figure [5] also includes the SLSN-II from the light curve 
sample. It is interesting to note that if the peak magnitude 
of SN 2006gy is corrected for host absorption, then the com- 
bined distribution of SLSN-I and SLSN-II peaks sharply at 

M 21.8 with SN2006tf the only SLSN-II falling be- 

low0. It is unclear if this event is the low luminosity tail 
of the SLSNe-II population or the bright tail of the "nor- 
mal" Typelln supernova distribution. Indeed, it is unclear 
if two such distinct populations exist in nature, or if there 
is simply a publication bias against Typelln supernovae of 
moderately bright luminosities. Any such distinction may be 
borne out by the larger samples of normal supernovae and 
SLSNe currently being gathered by searches such as PTF, 
Pan-STARRS, the Catalina Real-Time T ransient Survey 
ilDrake et all 120091 ). and La Silla-QUEST (|Hadiivska et all 
l2012t ). 

The discovery of such rare events, however, may be 
hindered by contamination from other transient phenom- 
ena. The ROTSE-IIIb sample includes about ten times more 
SNela than SLSNe-like events even though the SLSNe are 



2 SN 2006tf was included in the sample of lGal-Yam] (|2012| 1 based 
on its peak above —21 mag as suggested by preliminary ROTSE- 
IIIb estimates, but our revised analysis and dimming to match 
the filtered photometry places it below this cutoff 
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observable over a much larger volume because the ratio of 
the rates is even larger. There is also the worry that other 
rare transient phenomena, such as tidal disruption flares or 
an outburst from an AGN with a supernova-like light curve, 
could contaminate the sample. We potentially have one such 
source in our ROTSE-IIIb sample: Dougie. If Dougie proves 
to be a hydrogen poor supernova, our SLSN-I rate would 
roughly double. 

The SLSN-II p opulation may sim ilarly be polluted 
with AGN activity. iDrake et al.l l|201lf ) favor classifying 
CSS100217 as a superluminous supernova even though it 
radiated far more light than other SLSN published so far. 
CSS100217 does show features present in Typelln super- 
nova, but these are not uniq ue to supernovae and can also 
mark an AGN o utburst fsee lColgate fe Cameronl 1 19631 and 
lFilippenkolll989l for a discussion of the similarities between 
AGN and SNelln spectra). Since the host of CSS100217 is 
known to harbor an AGN and this is precisely coincident 
with the transient, this co uld simply b e a de monstration of 
rare AGN related activity. IDrake et al.l (|201lr ) note that am- 
plitude of CSS100217's outburst stands distinct from typical 
AGN, but CSS100217 was specifically selected for this un- 
usual brightening while the comparison sample was drawn 
from the SDSS DR3 spectroscopic sample (i.e. not selected 
based on variability). The relatively large effective volume- 
time for CSS100217-like events listed in Table indicates 
that such events are either very rare (least they dominate 
flux limited surveys, which is not the case), or they must be 
selected against due to, for example, the presence of a co- 
incident AGN. We know of no reason to exclude AGN from 
having supernova-like outbursts, so it will be important to 
measure the rate of such events to determine their potential 
contamination to supernovae samples. 

Like CSS100217, SN2006gy was located near the core 
of its host (but clearly offset in this case). Studies of our 
own galaxy and M31 have uncovered a significant popula- 
tion of massive star s located very close to the central nuclei 
(|Genzel et al.l |2000| ; lLauer et al.ll2012T V We speculate that 
such environments may on occ asion harbor dense clouds o f 
gas accreted from the disk (cf. lLevin fc Beloborodovll2003h . 
recent galaxy mergers, or fed by filaments, and these spe- 
cial environments may yet give birth to unusual supernovae. 
These events may be heavily extincted and thus difficult 
to detect (e.g. SN 2006gy suffers from about 1.25 mag of ab- 
sorption in the R-band). This may add to the confusion with 
AGN and make the SLSN-II rate difficult to measure. 

The rates presented here may be helpful in deter- 
mining the origin of SLSNe. For exa mple, using the to- 
tal star formation rate at z=0.16 (iHopk ins fc Beacoml 
20061') and a modified Salpeter IMF |Baldrv fc Glazebrook 
20031 ), we can compare our total SLSN rate to the ex- 
pected production rate of stars with 56Mq < Mi n i t i a i < 
84Mq, which is the mass range under which zero 
metallicity stars are exp ected to enter the puls a tional 
pair instability regime (|Chatzopoulos fc Wheelerl |2012|; 
lYoon. Dierks fc Langerll2012l ). We find that our SLSN rate 
could be explained by a small fraction of the stars born 
in this mass range (about 1 to 3% at the la confidence 
level). To put this another way, we can say at the > 10cr 
confidence level that at z ~ 0.2, not all stars in this mass 
range produce supernovae brighter than M ~ —21 mag; if 
some of these stars do produce SLSNe, others must expe- 



rience alternate fates. These rates also define a basis for 
comparison to future, higher redshift studies, which may 
help determine if the IMF varies over the distances that will 
be probed by forthcoming surveys (i.e. z ~ 4 with Subaru 
Hyper-SuprimeCam) . 
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